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Abstract 



a 

c/3 I Numerical investigations have been made of responses of a Hodgkin-Huxley (HH) neuron 

to spike-train inputs whose interspike interval (ISI) is modulated by deterministic, semi- 
deterministic (chaotic) and stochastic signals. As deterministic one, we adopt inputs with 
the time-independent ISI and with time-dependent ISI modulated by sinusoidal signal. The 
Rossler and Lorentz models are adopted for chaotic modulations of ISI. Stochastic ISI inputs 
with the Gamma distribution are employed. It is shown that distribution of output ISI data 



O ' depends not only on the mean of ISIs of spike-train inputs but also on their fluctuations. 

^ ■ The distinction of responses to the three kinds of inputs can be made by return maps of 

input and output ISIs, but not by their histograms. The relation between the variations of 
►>. ! input and output ISIs is shown to be different from that of the integrate and fire (IF) model 

O I because of the refractory period in the HH neuron. 
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I. INTRODUCTION 

Neurons in our brain are known to be responsible for encoding the characteristics of 
stimuh into a form for further processing by other neurons. During the last decades the 
anatomical, physiological and theoretical studies on neurons have been extensively made. 
Despite these efforts, the code used for encoding and decoding in neurons is not clarified at 
the moment |l[]. It is commonly believed that the firing rate reflects the strength of the inputs 
which trigger the action potentials of neurons. Indeed, the firing activities of motor and 
sensory neurons vary in response to the applied stimuli. It is not known, however, whether 
the information is carried through the mean firing rate (rate encoding) or through the 
details of sequences of the temporarily encoded interspike interval (ISI) (temporal encoding), 
which is currently controversial |^- [Q. In the last few years, experimental evidences have 
accumulated, indicating that many biological systems use the temporal coding. Human 
visual systems, for example, have shown to classify patterns within 150 msec in spite of the 
fact that at least ten synaptic stages are involved from retina to the temporal brain [Q . The 
similar speed of visual processing has been reported for macaque monkeys 0]. Because the 
firing frequency of neurons involved is less than 100 Hz, each neuron can contribute at most 
one or two spikes to such computations; there is not sufficient time to sample firing rates. 

In recent years, much studies on the encoding of the spike trains by neurons have been 
made by using the integrate and fire (IF) model 0], which is one of the simplest, dynamical 
models of neurons |§]. The IF neutron is silent without the external, input current, I\. 
When Jj, exceeds the critical value, /jc, the IF neuron shows the self-excited oscillations, 
whose frequency, /o, depends on the magnitudes of Jj. It is shown that /o continuously 
vanishes when Jj is decreased and approaches to lie. This behavior of the continuous /o — h 
dependence is different from the discontinuous one at lie in the more realistic Hodgkin-Huxley 
(HH) neurons 0; the IF and HH neurons are classified as the type I and type II, respectively 
|T0| . Furthermore, the IF neuron has the disadvantages of the artificial reset of the action 



potential and the lack of the refractory time. Although it has been widely employed for the 
study of neural networks, the IF model is too crude to discuss the activities of real neurons. 
The HH model, which well describes the spiking behavior and refractory properties of 
real neurons, is expressed based on non-linear conductances of Na and K ion channels 0. 
Since the HH model was proposed, its property has been intensively investigated 
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The behavior of self-excited oscillations of the HH neuron with the applied current has 
much variety than that of the IF model. It is shown that the oscillation of the HH neuron 
may become chaotic when the sinusoidal Ji is applied with proper choices of magnitude and 



frequency [|^ fl^. Such chaotic oscillations are experimentally observed in squid giant 
axons 1 14 1 |15| and Onchidium neurons ||18|| . 

The HH model was originally proposed to account for the property of squid giant axons 
and it has been generalized with modifications of ion conductances |jl9|. The HH-type models 



have been widely adopted for a study on activities of transducer neurons such as motor 
and thalamus relay neurons, which transform the amplitude-modulated input to spike-train 
outputs. In this paper, we pay our attention to data-processing neurons which receive and 
emit the spike-train pulses. Assuming that the data-processing neuron may be essentially 
described by the ion-conductance mechanism of the HH model, we investigate its input- 
output response in order to get some insight to the following questions: 



(1) How the output ISIs depend on the input ISIs? Does the average rate of the output ISI 
depend only on the average of the input ISIs? 

(2) How do neurons distinguish the different types of deterministic, chaotic and stochastic 
inputs? How different is the response to different types of spike-train inputs? 

Our paper is organized as follows: In the next sec. II, we mention a simple neuron model 
adopted for our numerical calculation. In sec. HI, we investigate the response of our system 
to deterministic inputs with time-independent ISI (§3.1) and time-dependent ISIs modulated 
by sinusoidal signal (§3.2). Input and output ISIs are studied by their histograms and return 
maps; the former shows the distributions and the latter the time correlation of ISI data. In 
sec. IV, chaotic inputs generated by Rossler (§4.1) and Lorentz model (§4.2) are discussed. 
Stochastic inputs with the Gamma distribution are treated in sec.V. The final section VI is 
devoted to conclusion and discussion 

II. ADOPTED MODEL 

We adopt a simple system consisting of a neuron and a synapse. The neuron is assumed 
to be described by the HH model and the synapse by the alpha function (Eq.(16)). We 
will investigate the response of our neuron when spike-train inputs are applied through the 
synapse. 

The HH model is described by the non-linear coupled differential equations for the four 
variables, V for the membrane potential, and m, h and n for the gating variables of Na and 
K channels, and it is given by 

CdV/dt = -g^,m^h{V - ^Na) - 9Kn\V - Vk) - 9l{V - Vi) + h, (1) 

dm/dt = — (om + bm) fn + a^, (2) 

dh/dt = -{ah + bh) h + ah, (3) 

dn/dt = -{an + bn)n + an, (4) 

where 

a^ = 0.1(l--f-40)/[l-e-(^+^°)/i°], (5) 

6^ = 4e-(^+65)/i«, (6) 

an = 0.01 {V + 55)/[l - e-(^+55)/i0|^ ^j^ 

6„ = 0.125 e-(^+6^)/^°, (8) 



6„ = l/[l + e-(^+=^^)/i°]. (10) 

Here the reversal potentials of Na, K channels and leakage are V^^. = 50 mV, Vk = ~77 
mV and l^ = —54.5 mV; the maximum values of corresponding conductivities are (^Na = 
120 mS/cm , g^ = 36 mS/cm and gi^ = 0.3 mS/cm ; the capacity of the membrane is 
C = 1 /iF/cm , details of the HH model being presented in Refs. p [^ . 



The external, input current, Ji, is taken to consist of two terms; 

/i = /s + /p. (11) 

where Is expresses the static dc current and Jp denotes the pulse current induced by the 
spike-train input whose explicit form will be discussed shortly (Eq.(15)). 
We consider the delta- function- type spike-train input expressed by 

Ht) = VaY.^it-tin). (12) 

n 

The firing time ti„ for arbitrary n is assumed to be recurrently defined by 

tin+l ^ tin ~r Iin\tin), V "-'J 

til = 0, (14) 

where ISI of input spike, Ti„, is generally a function of a given time ti„. In this study, we 
take Tin to be constant ISI, and time-dependent ISI modulated by sinusoidal, chaotic and 
stochastic signals. 

The spike train given by Eq.(12) is assumed to be injected through the synapse, yielding 
the current Jp given by 

Ipit) = gsyn J2 ^i^" ^i") (^'^ ~ ^syn)- (15) 

n 

Here f^syn and Vgyn are the conductivity and potential of synapse, and the alpha function, 
a(t), is defined by [^ 

a{t) = it/r) e-''^ e(t), (16) 

where r is the time constant relevant to the synapse conduction and G(t) is the Heaviside 
step function. When the ISI is very large compared with r, Eqs.(15) and (16) yield pulse 
currents with the maximum value of J™*^ = e~^ gsyn (K — V^yn) at t = ti„ + r and with the 
half-width of 2.45 r. We assume V^ = 30 mV (the typical value of the maximum membrane 
potential), Kjyn = —50 mV and r = 2 msec, and treat g^yn as a parameter. 

When the membrane potential V oscillates, it yields the spike-train output, which may 
be expressed by 

Uo{t) = VaY. S{t-tom), (17) 

m 

in a way similar to Eq.(12), and the output ISI is defined by 



We will investigate how Tom depends on the various types of Ti„. 

Differential equations given by Eqs.(l)-(10) including the external current given by 
Eqs.(ll)-(16) are solved by the forth-order Runge-Kutta method for 20 sec with the in- 
tegration time step of 0.01 msec. We discard results of initial ten thousand steps to get 
stable solutions. If ISI of spike-train input or output is about 10 msec, the size of its sample 
is about 2000. Although this figure is not sufficiently large for statistics of ISI data, we hope 
an essential ingredient will be clarified in our numerical investigation. 



III. DETERMINISTIC INPUTS 

3.1 Time-Independent ISI 

A. Pacemaker Neurons 

Let us first consider the HH neuron without the spike-train input (Jp = 0). The HH 
neuron is reported to be silent for Jg = 0, and to show the self-excited oscillation when Ig 
exceeds the critical value of lie = 6.3 /xA/cm , above which Ton decreases gradually as Ig is 
increased. The dashed curve in Fig. 1(d) expresses an example of the self-excited oscillation 
with the period of Ton = 10.75 msec for Is = 25 /xA/cm and Jp = 0. 

Now we apply the spike-train input to this self-excited (pacemaker) neuron. The input 
is given by Eqs.(10)-(16) with the time-independent ISI of Ti„(tin) = Ti = 15 msec and 
fi'syn = 0.5 mS/cm . This spike-train input, Ui, shown in Fig.l(a) yields the pulse current, 
Jp, shown in Fig. 1(c), by which the membrane potential, V, oscillates as depicted by the 
solid curve in Fig. 1(d). We plot in Fig. 1(b) the time sequence of the spike-train output, 
Uo, which should be compared with the input Ui. The pulse current, Jp, has the maximum 
value of J™*^^ = 14.8 yuA/cm at t = ti„ + r msec. We notice that the oscillation in V is 
rather different from that shown by the dashed curve for Jp = 0. Figure 2 expresses the 
histogram of the output ISI, showing that {To} distributes continuously between 8.36 to 
11.62 msec. The mean and root-mean-square (rms) values of the output ISI are fio = 10.43 
and (Jo = 1.12 msec, respectively. This oscillation is chaotic as was pointed out for the HH 
neuron receiving sinusoidal inputs [|16| Ijl^; the pulse current Jp shown in Fig. 1(c) is not so 



different from the sinusoidal one in a crude sense. The chaotic behavior is clearly seen in 
Fig. 3, which depicts return maps of input and output ISIs. 

When the ISI value of spike-train input, Tj (= /ii), is changed, we get an interesting 
behavior in /Iq as shown in Fig. 4, where the solid (dashed) curve expresses /Xq (o'o), and 
filled circles express the distribution of {To„} for a given /ii. We note that for /ii = 9 — 11 
msec, the period of the oscillation is forced to be the same; Tom = Ti, leading to the ratio: 
k = /io//ii = 1. When /ii = 5 msec, we get /io = 10 msec and then k = 2. On the contrary, 
for /ii = 20 msec, we get the two values of Ton = 9 and 11 msec, and /io = 10, the average 
period of the output being a half of the input {k = 1/2). This is also the case for /ii = 21 
and 22 msec. In the other cases noticed above, the ISI of output distributes between about 
8.5 - 11.5 msec. We should note that irrespective of /ij, output ISI is always about 10 msec, 
which is nearly equal to To = 10.75 msec, ISI for Jg = 25 and Jp = /iA/cm . 



B. Silent Neurons 

Next consider tlie silent neuron with Jg = 0, for which the oscillation of the membrane 
potential is induced by applied spike-train inputs. Figures 5(a)-(d) show the calculated 
result in which the spike-train input is given by Ti„ = yUi = 10 msec and Qsyn = 0.5 mS/cm 
without static currents (Is = 0). The applied spike-train inputs shown in Fig.5(a) create 
the pulse current with the peaks of J™*^^ = 15.3 /xA/cm as shown in Fig. 5(c). The induced 
oscillation of the membrane potential, V, in Fig. 5(d) is phase locked with the ratio of 4 : 3, 
oscillating with a long cycle of 40.00 msec (=11.25 + 12.36 + 16.39) = 4 /ij, where 11.25, 
12.36 and 16.39 are the values of output ISIs. The return map of output ISIs is plotted in 
Fig.6(a). 

Figure 7 shows /Xq and cTq as a function of /ii. We notice that /Iq agrees with /i; {k = 1) for 
Hi greater than 12 msec, where the HH neuron behaves as a simple transmitter with a delay 
of about 2.0 msec. This is in strong contrast with the behavior of the pacemaker neuron 
discussed in the preceding sub-section (Fig.4). On the other hand, for /ii less than 11 msec, 
the behavior of output ISI is rather complicated. It is easy to see that /c = 2 for /Xj = 6, 7 
and 8 msec, and that fc = 3 for /ii = 4 msec. For /ii = 9, we get T^m = 12.06 and 14.96 msec, 
leading to a longer period of 3/ij = 27.00 ( = 12.06 + 14.96) msec. For /ii = 5 msec, we get 
Tom = 10.94 and 14.06 msec, which leads to a long period of 5 /ii = 25 msec, its return map 
being shown in Fig.6(b). Surprisingly, a much longer period of 13 /Xi is realized for /ii = 11 
msec. The rms value of o"o has an appreciable value only around /ii = 10 msec. 

We have repeated our calculation by changing the value of (/syn- The calculated ratio, 
k = /io//^i, is shown as functions of gsyn and /ii in Fig. 8, where only the integer values of 
k are shown by circles {k = 1), squares(/i; = 2), triangles {k = 3) and diamond {k = 4). 
Note that non-integer values of k exist between the integer values; for example, k = 4/3 
for gsyn = 0.5 mS/cm and /ii = 10 msec (Fig.7). We cannot obtain spike-train outputs 
for small synaptic couplings as expected. When /ii = 10 msec, we get the critical value of 
fl'syn = 0.11 mS/cm below which no outputs are available. This coupling yields the pulse 
current with the maximum value of J™^^ = 1.6 /iA/cm , which is much smaller than the 
critical dc current of lie = 6.3 /iA/cm for the self-excited oscillation with Jp = 0. We note 
that we get k = 1 for the large ISIs with fairly strong synaptic couplings. When we decrease 
/ii with keeping ^fgyn fixed, values of k become larger since the HH neuron cannot respond 
to inputs with the small ISI because of its refractory period. Figure 8 reminds us the result 
of Guttman, Feldman and Jakobson |jl5| who reported in their Table 1, the calculated k as 
functions of the magnitude. A, and the frequency, /i, when the sinusoidal input given by 
I[ = A sin(27r /i) + Jb is applied to squid giant axons with a bias current, Jb. Our result for 
/ii < 10 msec fairly agrees with that of Ref. [|T3[. However, the agreement between the two 
results is not good for /ii ^ 10 msec, where our input current with the pulse width of about 
2.45 r ~ 5 msec (Eqs.(15) and (16)) is quite different from the sinusoidal current adopted 



in Ref. 15 



As was shown in Fig.4, the pacemaker HH neuron emits the output ISI of Tom ~ 10 msec 
irrespective of the value of input ISI, and then it is considered to be inadequate as a data 
processor. Then, in the following sections, we will investigate only the silent HH neuron 
with a fixed value of gsyn = 0.5 mS/cm . 

3.2 ISI with Sinusoidal Modulation 



In this sub-section we discuss an application of the spike-train input whose ISI is mod- 
ulated by the sinusoidal signal given by 

Ti„(t) = rfo + rfisin(27rt/Tp), (19) 

where Tp is the period and do and di are coefficients adjusting fi^ and o"i. 

Figure 9(a) and (b) show the time course of input Ui and output f/o for d^ = 2di = 
20 msec with Tp = 100 msec. Because of the introduced sinusoidal modulation, ISIs at 
100 < t < 150 msec are, for example, larger that those at 150 < t < 200 msec. Figure 
9(c) depicts the return map of input ISIs, which has the egg-shape circle expected for the 
sinusoidal signal. On the other hand, the return map of output ISIs shown in Fig. 9(d) 
reveals the chaotic behavior. Results for d^ = 2di = 10 are plotted in Fig. 10 (a)- (d), which 
show that although the return map of input ISI has the egg-shape circle, that of output 
ISIs is distorted. The reason of this distortion is explained in Fig. 11, where solid histograms 
express input and output ISIs for do = 2di = 10 (/i = 8.68, a = 3.42) and dashed histograms 
those for do = 2di = 20 (/Xj = 17.54, ctj = 6.94), with Tp = 100 msec. (It is noted that 
we get fj,i < do because the histogram of the input ISI at Ti„ < do has larger magnitudes 
than that at Ti„ > do). In the case of do = 2di = 20 msec, the input and output ISIs 
distribute almost in the same region at 11 < Tom < 30 msec. On the contrary, in the case 
of do = 2di = 10 msec, the output ISIs distribute at 11.01 < Tom < 19.48 msec while input 
ISIs are at 5.00 < Ti„ < 14.96 msec; no output ISIs at Tom < H msec. This is due to 
the refractory period of the HH neuron and it is the origin of the distortion in the return 
map shown in Fig.lO(b). Defining the dimensionless coefficients of variations for input and 
output ISIs by 

CvA = crx/fi\, (A = i and o) (20) 

we get Cvo = 0.17 and 0.38 for do = 2di = 10 and 20 msec, respectively; note that c^i = 0.40 
for both inputs. 

Figure 12 shows /Iq and (Tq calculated by changing /ij with the fixed value of Cvi = 0.40. 
Solid and dashed curves denote fio and (Tq, respectively, and filled circles the distribution of 
{Tom} for a given /x;. We notice that there is no output ISIs with Tom less than about 10 
msec, which shows characteristic of the low-pass filter of the silent HH neuron . 



IV. CHAOTIC INPUTS 

4.1 Rossler Model 

In this section, we study the spike-train input whose ISI is modulated by chaotic signals. 
First we adopt the Rossler model, which is given by 

dx/dt = —y — z, (21) 

dy/dt = X + ay, (22) 

dz/dt = bx — c z + xz, (23) 



with a = 0.36, b = 0.4 and c = 4.5 |2T|. Since ISI has to be positive and the characteristic 
time scale in the Rossler model is different from that of the HH model, we adopt the variable 
x{t) which yields 

TUtin) = do + (t/i/10) x{pUn), (24) 

with the following two choices of parameters: 

do = di = 10 msec and p = 1/10 (case Rl), 
do = di = 20 msec and p = 1/20 (case R2). 

Figure 13(a) and (b) show the time course of input and output spike trains for the 
case Rl. The return map of input ISI depicted in Fig. 13(c) shows a shape characteristic 
for chaotic signals. On the other hand, the return map of output ISIs is rather strange 
with no traces at Tq < 10 msec. This is due to the low-pass filter behavior of the silent 
HH neuron, as shown by solid histograms in Fig. 14(a) and (b); output ISIs distribute at 
11.11 < Tom < 25.15 msec {fiom = 13.43, o"o = 2.44 msec) whereas input ISIs distribute at 
5.06 < Tin < 16.56 msec (/ii = 9.53, cri = 2.69 msec). 

Return maps for the case R2 are shown in Fig. 15, in which both return maps are almost 
the same. This is because input and output ISIs locate almost in the same region of 10 < 
Tin, Tom < 30 msec, as shown by dashed histograms in Fig. 14(a) and (b). 

Next we investigate the nature of the correlation in the ISI sequences. This is made 



by employing the surrogate data method applied to ISI data [23]. We adopt the shuffled 
surrogate as a simple method to get surrogate data. The distributions of ISIs of shuffled 
surrogate inputs are exactly the same as those of original ISI data although surrogate data 
have no correlation between successive ISI values. 

The time course of the membrane potentials for the surrogate data is ostensibly quite 
similar to that for the original chaotic input (not shown). The solid (dashed) histogram in 
Fig. 14(c) shows the distribution of output ISIs of surrogate data generated from the Rossler 
model for the case Rl (R2). The results of the surrogate data are similar to those for the 
corresponding original data. Return maps of the input and output ISI of the surrogate data, 
depicted in Fig. 15(c) and (d), show the characteristics of random signals. 

4.2 Lorentz Model 

The similar calculation is made with the use of the Lorentz model, which is given by 

dx/dt = d {y — x), (25) 

dy/dt = ex — y — xz, (26) 

dz/dt = -f z + xy, (27) 

with d = 10, e = 28 and / = 8/3 p2[. We employ the variable z(t), with which the input 
ISI is given by 

Tin(t„) = do + (rfi/25) iz{ptin)-25), (28) 



where do = di = 20 and p = 1/100. 

Figures 16(a)-(d) show return maps of ISI data of original chaotic data and its surrogate. 
Return maps of output ISIs for the chaotic and surrogate data shown in Fig. 16(b) and (d) 
have no traces at T^m > 10 msec because of the low-pass filter character of the HH neuron. 

Sauer ||2^, and Racicot and Longtin |^ studied the response of the IF model to the 



input whose amplitudes are modulated by chaotic signals. It was shown that when the mean 
firing rate is high, the relationship between input and output is high, which leads to the 
high nonlinear predicability. Our calculations have shown that when the mean firing rate 
is too high [i.e. input ISIs are too short such as Tin < 10 msec), the information is lost 
because the HH neuron behaves as the low-pass filter due to its refractory period, which is 
not included in the IF model. 



V. STOCHASTIC INPUTS 

The ISIs of spike-train input. Tin, in Eq.(13) are assumed to be independent random 
variables with the Gamma probability density function given by 

P{T) = s' T'-^ e-'^l r(r) (29) 

for which we get /ii = r/s, Oi = \frls and c„j = l/v^j T (r) being the gamma function. It 
is noted that in the limit of r = 1, Eq.(29) reduces to the Poisson distribution (cvi = 1) and 
that in the limits of r ^ cxd and s — > oo with keeping /ii = rjs fixed, Eq.(29) reduces to 
PiT^ = S(T — /ii), the constant ISI with fii = T and Cvi = 0. 

The spike-train input created by the Gamma-distribution generator is applied to our 
neural system. Calculations are performed by changing /ii by keeping the value of Cvi fixed. 
Note that because the number of our sample of input ISI is not sufficiently large, the obtained 
Cvi fluctuates around the intended values. Solid histograms in Fig. 17(a) and (b) show the 
result for Cvi = 0.40, /ii = 10 msec, Cvo = 0.25 and /Xo = 14.84 msec while dashed histograms 
for the result for Cvi = 0.40, /ii = 20 msec, Cvo = 0.36 and /io = 21.11 msec. 

Solid and dashed curves in Fig. 18(a) denote /Xq and cTq, respectively, for Cvi = 0.4. We 
note that as increasing /Xi, /io increases and approaches the dotted line expressing /Xq = /Xi. 
This is similar to the case of c^i = shown in Fig. 7, where /Xq = /ii at /Xj ~ 10 msec. On the 
contrary, the dependence of /Xo on /Xj for the case of Cvi = 1.0 shown in Fig. 18(b), is rather 
different from the cases of cVi = and 0.4. We get /io ~ (/ii + 10) msec at /io < 100 msec and 
it deviates from the dotted line showing /io = /ii. These calculations depicted in Figs. 18(a) 
and (b) clearly show that /io depends not only on /i^ but also on ai (c^i). Although both /io 
and (To increase as the value of /ii is increased, the increase in the latter is more significant 
than that in the former, which yields an increase in Cvo, as shown by the thin-solid curves 
in Figs. 18(a) and (b). We note that Cvo approaches the values of Cyi as increasing /Xi, the 
related discussion being given in Sec. VI. 

We have performed the calculation also using input ISI with random, uniform distribu- 
tion. Obtained results are similar to those for the Gamma distribution, as far as the adopted 
values of Cvi are the same (not shown). 



VI. CONCLUSION AND DISCUSSION 



We have investigated the responses of the HH neurons, by applying the various types of 
spike-train inputs whose ISI is modulated by deterministic, chaotic and stochastic signals. 
The obtained results are summarized as follows: 

(1) Output ISIs of the pacemaker HH neuron against the time- independent input ISI are 
always Tom ~ 10 msec irrespective of the value of the input ISI (Fig.4). 

(2) Output ISIs of the silent HH neurons for the constant ISI with T^ > 10 msec yield 
output ISI with Tom = Tin whereas for ISI with Ti„ < 10 msec, the HH neuron generally 
emits multiple numbers of output ISIs (Fig. 7). 

(3) For the input ISI modulated by sinusoidal, chaotic and stochastic signals, the silent HH 
neuron behaves as a low-pass filter because of its refractory period, yielding output ISI with 
Tom > 10 msec. 

(4) Output ISIs generally depend not only on the mean of the input ISI but also on their 
fluctuations: the HH neuron is not a simple integrator. 

(5) The analysis on the histograms of input and output ISIs cannot distinguish the responses 
to the deterministic, chaotic and stochastic signals. 

(6) The distinction can be made by an analysis of the time correlation of the ISI data, for 
example, by plotting their return maps. 

Softky and Koch [EBl have reported a large coefficient of variation (cvo = 0.5 ~ 1.0) 



for spike trains of non-bursting cortical neurons in visual VI and MT of monkeys in strong 



contrast with a small c^o (= 0.05 ~ 0.1) in motor neurons |2^. In order to explain the large 
Cvo, several hypotheses have been proposed; a balance between excitatory and inhibitory 
inputs [^, the high physiological gain in the /o — /i plot ||2^, correlation fluctuations in 



recurrent networks [0, and the active dendrite conductance |^T|]. By using the IF model 



Feng and Brown |^2|] have shown that there are three kinds of behaviors of Cvo depending 
on the distribution of input ISIs; 

(a) Cvo tends to decrease for the Gaussian, uniform or truncated distribution of ISIs, 

(b) Cvo remains constant for the exponentially distributed ISIs, and 

(c) Cvo diverges to infinity when ISIs follow the Pareto distribution which has a slow- 
decreasing tail of T~" (a > 0) at large T. 



The case (a) was previously discussed by Marsalek, Koch and Maunsell |^ . 

Figure 19 shows the dependence of Cvo on Cvi for various types of our input ISIs having 
been reported in previous sections. Inverted triangles denote the results for the constant ISIs 
(§3. IB), open marks the results for input ISIs with sinusoidal modulation (§ IV), and filled 
circles, triangles and squares results of the stochastic modulation (§ V). Our calculations 
show the followings: 
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(i) The constant ISI with a vanishing Cvi yields the finite Cvo (< 0.2), i.e. Cvo ~ c 

(ii) the finite-width distribution of ISIs with the sinusoidal modulation leads to Cvo ~ Cvi, 

and 

(iii) the exponential. Gamma distribution of input ISIs yield the result which is ostensibly 

similar to that in the item (ii). 

Although the item (ii) is in agreement with the result for the above-mentioned case (a) 

discussed in Refs. |^ and |^ , the items (i) and (iii) disagree with the results of the cases 

of (a) and (b), respectively. This difference is expected to arise from the fact that the 

response of the HH neuron with the refractory period is different from that of the IF neuron 



without it [ p!0[] . 
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Finally we want to discuss the transient response of the HH neuron to the cluster of 
spike-train inputs. Figure 20(a), (b), and (c) show the results for Ti„ = 5, 10 and 20 msec, 
respectively. In Fig. 20(c), for example, the upper (lower) panel of CI, C2, C3 and C4 
express the time courses of input (output) spike trains for inputs of two, three, four and five 
impulses, respectively, with Tin = 20 msec. In this case the ISI of output pulses is the same 
as that of input pulses, Tom = T^, because the HH neuron behaves as a linear transmitter 
for inputs with ISI of Tm ~ 10 msec. On the contrary, its behaviors of output ISI data for 
inputs with T[n = 5 and 10 msec are much complicated. Fig. 20 (a) shows that output ISIs 
for Tin = 5 msec are 11.39 and 11.87 msec, which should be compared with T^m = 10.94 and 
14.06 msec for the sequence of the spike trains with constant ISI of 5 msec discussed in sec. 
3. IB. For the case of Ti„ = 10 msec shown in Fig. 20(b), we get Tom = 11.44, 11.80 and 17.11 
msec whereas the sequence of the constant ISI of 10 msec leads to Tom = 11.25, 12.36 and 
16.39 msec (Fig.5). It is noted that both inputs with three (A2) and four impulses (A3) yield 
the same output of two impulses. Similarly, inputs with three (B2) and four impulses (B3) 
lead to outputs with three impulses. We should note in all the cases shown in Fig. 20 (a)- (c) 
that the first output pulse is rather quickly emitted with a delay of 2.1 msec after the first 
input pulse of clusters is applied to the HH neuron. This fast transient response may be 
relevant to a quick passage of information reported by Thorpe, Eize and Marlot and by 
Rolls and Tovee flBI. 
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FIGURES 



FIG. 1. Responses of the pacemaker HH neuron to the time-independent input ISI (/p with 
Ti = 15 msec and Is = 25^A/cm ; time sequences of (a) input [/;, (b) output Uq, (c) pulse current 
/p, and (d) membrane potential V, dashed curves in (b) and (d) denoting the result with /p = 0. 



FIG. 2. Histograms of (a) the time-independent input ISI (Tj = 15 msec) and (b) output ISI 
of the pacemaker neuron (see Fig. 1). 

FIG. 3. Return maps of (a) time-independent input ISI with Tj = 15 msec and (b) output ISI 
of the pacemaker neuron (Fig.l). 



FIG. 4. Mean (^Uq, solid curve) and rms (do, dashed curve) values of output ISI of pacemaker 
neurons (/g = 25^ A/cm ) against the mean value (/ii) of time-independent input ISI. Filled circles 
denote the distribution of output ISIs for a given ^Uj, dotted curves denoting k = /Uo//^i. 



FIG. 5. Response of the silent HH neuron to the time-independent input ISI (/p with Tj = 10 
msec and /g = 0); time sequences of (a) the input Ui, (b) output Uq, (c) pulse current /p, and (d) 
membrane potential V. 

FIG. 6. Return maps of output ISIs for the time-independent input with (a) T[ = 10 and (b) 
5 msec (see Fig. 5). 

FIG. 7. Mean (/io, solid curve) and rms (o"o, dashed curve) values of output ISI of silent neurons 
against the mean value (/ii) of time-independent input ISI. Filled circles denote the distribution of 
output ISIs for a given /ii, dotted curves denoting k = /io//Ui- 

FIG. 8. The calculated ratio, k = /Uo//ii, as functions of /ii and gsyn for the time-independent 
ISI input to silent neurons: only integer fc's are shown by circles (fe = 1), squares (A; = 2), triangles 
(A; = 3) and diamond (A; = 4), the cross denoting no outputs. 

FIG. 9. Time courses of (a) input Ui and (b) output Uq, and return maps of (c) input and (d) 
output ISIs for the sinusoidal modulation for do = 2di = 20 msec (Eq.l9). 

FIG. 10. Time courses of (a) input Ui and (b) output Uq, and return maps of (c) input and 
(d) output ISIs for the sinusoidal modulation for do = 2di = 10 msec (Eq.l9). 

FIG. 11. Histograms of (a) input ISI with sinusoidal modulation and (b) output ISI. Solid 
(dashed) curves are for do = 2di = 10 (20) msec. 
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FIG. 12. Mean (po, solid curve) and rms (do, dashed curve) values of output ISI against the 
mean value (//i) of input ISI with sinusoidal modulation (cvi = 0.40). Filled circles denote the 
distribution of output ISIs for a given /ii, dotted curves denoting k = fio/fii. 

FIG. 13. Time courses of (a) input Ui and (b) output Uo, and return maps of (c) input and 
(d) output ISIs for chaotic inputs generated by the Rossler model (case Rl). 

FIG. 14. Histograms of (a) input ISI and (b) output ISI for the chaotic input generated by the 
Rossler model, and (c) output ISI for its surrogate, solid (dashed) curves being for case Rl (R2). 

FIG. 15. Return maps of (a) input ISI and (b) output ISI for the chaotic input generated by 
the Rossler model (case R2); (c) and (d) are corresponding return maps of its surrogate data. 

FIG. 16. Return maps of (a) input ISI and (b) output ISI for the chaotic input generated by 
the Lorentz model; (c) and (d) are corresponding return maps of its surrogate data. 

FIG. 17. Histograms of (a) input ISI and (b) output ISI for spike-train inputs with the Gamma 
distribution; solid (dashed) curves for input ISIs of fii = 10 (20) msec with Cvi = 0.40. 

FIG. 18. Mean (^o) solid curves), rms ((Tq, dashed curves) and c^o (thin solid curves) of output 
ISIs against the mean value {p,i) of input ISI for Gamma distribution with (a) Cvi = 0.4 and (b) Cvi 
= 1.0. The dotted curves denoting fio = Hi are plotted for a guide of the eye 

FIG. 19. Cvo against c^[ for inputs with time-independent ISI (inverted triangles, Cvi = 0), 
with sinusoidal modulation (open diamonds, squares and triangles for Cvi = 0.22, 0.40 and 0.82, 
respectively), and with random Gamma distribution (filled squares, triangles and inverted triangles 
for Cvi = 0.40, 0.75 and 1.03, respectively). The dotted curve expressing Cvo = Cvi is plotted for a 
guide of the eye. 

FIG. 20. Time courses of spike-train inputs and outputs; input ISIs are (a) Tj = 5, (b) 10 and 
(c) 20 msec, and upper (lower) panel of each figure shows inputs (outputs). 
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